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Abstract 

o 

^vj , The exact results for thermodynamical properties of one-di- 

f"^ ■ mensional spin-i isotropic XY model with Dzyaloshinskii-Mo- 

£ — ' riya interaction in random lorentzian transverse field are obtained. 

0^ ■ This permits to discuss some approximate methods of disordered 

spin systems theory. The approximate scheme of examining the 
thermodynamics of one-dimensional spin-^ XXZ Heisenberg mo- 
del with Dzyaloshinskii-Moriya interaction in random lorentzian 
field is suggested. 
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Shortly after famous paper by E.Lieb, T.Schultz and D.Mattis El, who 
making use of Jordan- Wigner transformation reformulated the Hamilto- 
nian of spin s = h XY chain in terms of non-interacting fermions and 
obtained a series of rigorous results, the random versions of such models 
attract much attention. Some exact results were derived by H.Nishimori 
[El for isotropic XY chain in random lorentzian transverse field. It be- 
came possible because of the fact that after fermionization of such model 
one comes to a system of electrons on lattice that may transfer from site 
to site with the random (lorentzian) energy at sites. The average one- 
particle Green's function for such model was obtained first by P.Lloyd 
§ (see also |,1, §§). 

What follows is based on the notation that similarly to M one can 
consider the case of more complicate interspin interaction including to it 
the so-called Dzyaloshinskii-Moriya interaction. It was introduced phe- 
nomenologically by I.E.Dzyaloshinskii B and then derived by T.Moriya 
[Qj and it is widely used as one of microscopical mechanism (together 
with ANNNI model) of appearance of incommensurate phase in crys- 
tals @. V.M.Kontorovich and V.M.Tsukernik noted Q that taking 
into account of such interaction in s = 5 XY chain does not destroy 
the consideration proposed in Ref. jjj since after Jordan- Wigner trans- 
formation, as in the previous case, one comes to a quadratic in Fermi 
operators form. In Ref. |0] the problem about the possibility of ap- 
pearance of spiral spin structure in such a chain was examined; for this 
purpose the pair equal-time spin correlation functions were evaluated. 
Other papers dealing with statistical mecanics of s — h XY chains with 
Dzyaloshinskii-Moriya interaction EJ, [l| Q {jj| [H| ^ [j]| || like the 
Ref. Ill) are devoted to perfect (non-random) versions of the model. 

The present paper contains some exact results of statistical mechanics 
of ID s = \ isotropic XY model with Dzyaloshinskii-Moriya interaction 
in random lorentzian field. The paper is organized as follows. Fermion- 
ization and Lloyd's problem are considered in Section 1. Here the av- 
erage one-fermion Green's functions, the average spectral density and 
the average fermion correlation functions are derived. Section 2 contains 
calculations of thermodynamical properties of the model in question. 
A discussion of the estimation of static spin correlations is also given 
in it. The comparison of exact results with the ones obtained within 
different approximate approaches (Bose commutation rules approxima- 
tion, Tyablikov-like approximation, coherent potential approximation) 
are performed in Section 3. The developed in Sections 1 and 2 scheme 
may be used for the approximate study of ID s = \ XXZ Heisen- 
berg model with Dzyaloshinskii-Moriya interaction in random lorentzian 



field. This possibility is discussed in Section 4. Conclusions are given 
in Section 5. Briefly the results of the present paper were reported in 



1 Fermionization, Lloyd's problem, average 
one-fermion Green's functions, aver age spec- 
tral density, average fermion correlation 
functions 

A chain of N spins s — | with interaction between nearest neighbours, 
that are in transverse fields with random component distributed accord- 
ing to lorentzian law is considered. The Hamiltonian of the model has 
the form 

N N-l 

H = £(fio + n^ + J E ( s * s Ui + s M+i) 

+ D E K-^i - * v i s U) 

i=i 

N 



3 = 1 V J 



^fJ + iD,_ J-iD _ , \ 

+ E ("a— s Pj + i + -^~ s i s t+i ) • ( L1 ) 

After Jordan- Wigner transformation 

ci = sT, Cj = {-2sl){-2 S z 2 )...{-2 S z J _ l )s- , j = 2, ...,7V, 

one comes to the following quadratic in Fermi operators Hamiltonian 

N 



c+ = at, cf = (-2af)(-2al)...(-2«|_ 1 )s+, j = 2, ...,7V, 

{c+, c,-} = % {c+, c+} = 0, {c 4 , c,-} - (1.2) 



H = ^(n o + n,-)(ct Ci -0 



i=i 

AT-1 



j=i v 



that can be treated like in Ref . H . 

Let's introduce the following retarded and advanced temperature 
two-times Green's functions |E3] 

Gl m (t) = TiO(±t) < {c n (t), c+ (0)} > . (1.4) 

The goal of further consideration is to find the average Green's functions 
Gnm (t) , where the average means 

/+oo /*+oo 

- oo •/ — oo 

w 1 r 

^^'••■)=n-^Tr^ (L5) 

(i.e. OjS are independently distributed according to the lorentzian prob- 
ability distribution density centred at Qj = with the width T). 
It is easy to get the equation of motion for (1.4), namely 

ijGlmit) = S(t)S nm + (fio + n»)G? m (t) 
+^^ +M „(i) + J -^-G*_ 1>m {t). (1.6) 

Using spectral representation of Green's functions (1.4) 

1 f + °° 
GZm(t) = 7T du>exp(-iwt)GZ m (u>), 

27T J-oo 

/+oo 
dte X p(iu>t)G* m (t) (1.7) 

-OO 

the equations (1.6) can be rewritten in the form 

wGL(«) = t + (no + n n )G£„» 

+ — 2~ G £n,mM + -2-d,m(«). (1-8) 

or in the form that is initial for locator expansion 
G+ m (w) = , . + , . • (1.9) 

LJ—{llO + il n ) LJ—{llo + il n ) 



Formal solution as a series with respect to intersite interaction reads 

On 777. 



G£ m (w) = 



w - (O + 0„) 



J + iD 


1 


On+l,m 


2 w- 
J -iD 


1 


- (Oq + ^n+l) 

On— 1,777, 



W - (fi + fi„) W - (fio + ^n-l) 



(1.10) 



It is easy to average GJ m (w ± ie) presented as (1.10). Really, one 



should perform the averaging of [ „ ±ie _ {ih+il:i) 
integral 



p, i.e. to calculate the 



+oo 
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Jo2 



ft? + r 2 
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do. 



uj±ie — (fio + fij) 

r 



J (% + ir)(fij - ir)[w±ie - (Q + ttj)] k i 



I--(l-ll) 



For retarded (advanced) Green's function the integrand has one pole 
in the lower (upper) half of the complex plane Qj and two poles in the 
upper (lower) half of the complex plane Qj . Thus expanding the contour 
of integration in the lower (upper) half of the complex plane and using 
the residuum theory one ends up with 






uj±ie — (Qq + Qj 



1 






w±ie — (fio T *r) 



(1.12) 



that is the famous property of lorentzian distribution. In result the 
average series (1.10) can be summed up with the result 



G% m {u) ±ie) = 



2 ^n+l,m( w ±ie 



u>±ie — (Oo T *r) 



uj±ie — (51q T *r) 



(1.13) 



Since the average Green's functions are translationally invariant the 
equations (1.13) can be solved with the help of transformation 

1 



G% m (L0 ±ie) = -j- ^exp[i(n - to)k]G;F(w ± ie), (1.14) 



K = 2™,„ = -f , -f + 1, ..., f - 1 for even N or n = -»f±, -»jl 



1, ..., ^2"^ for odd N; the solution of algebraic equation for Gj(w ± ze) 
reads 



G^(co±ie) 



1 



w - [fi + V J 2 + £> 2 cos (k + p)]±i(e + T) ' 



costp : 



J 



VJ 2 + D 2 



, sirup 



D 



V J 2 + z? 2 



(1.15) 



The average Green's functions in site representation can be found in re- 
sult of summation over k in (1.14) that in thermodynamical limit reduces 
to the following integration 



1 f 

GZ m (u) ± it) = — / dn cxp \i(n — rn)K] 

27r J -IT 



u-{Q a + VJ2 + g2 {exp [»(« + p)] + cxp [-»(« + p)]} ± i(e + T)) 



.(1.16) 



Setting z — exp [i(n + </>)] one comes to the contour integral on unit 
circle in complex plane z 



Gt m (u> ± it 

x (h dz 



cxp [iip(n — m)] 

2~7ri 



&E±m z 2 _ [ W _ fio± , (e + r)]z + ^ZH2E 



(1.17) 



The denominator in the integrand in (1.17) should be presented as a 
product a(z — Z\)(z — z-z) where 



z\ = 



zi = 



Vj 


2 + D 2 
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- flo±i(e 


+ r) 




VJ 2 + D 2 


LU - 


- flo±i(e 


+ r) 



cu-n ±i(e + T) 



VJ 2 + D 2 



-1, 



V J 2 + D 2 



LU-tt ±i(e + T) 



1 2 



V J 2 + D 2 



-1, (1.18) 



and since z±Z2 — 1 z\[zi) is over (in) the unit circle. For n > m the 
integration yields 



exp [i 



ru-Oo±i(£+r) _ /r ^-Oo±i(e+r) 12 7l r 



V J 2 + D 2 



^-n ±i(e+r) 12 _ i 



41.19) 



whereas for n < 



f U — i1 ±i(e+r) /r ^-O ±t(e+r) i2 i 1 m-n 

exp [itp(n - m)\ \ ^p+m V [ yp+p? j 1 J ^ 2Q , 

V J 2 + D 2 / r o;-no±i(e+r) 1 2 \ 

Combing (1.19), (1.20) one finally gets for (1.17) 



Gtm{io ±ie) = 



f uj-n ±i(e+r) _ /r u>~n ±i(e+r) i 2 ~~p \n-m\ 

exp [iip(n - m)\ \ ^jr+W V [ V.P+D? J > ^ ^ 

V [ VJ 2 +D 2 J 

The obtained result (1.21) gives the average elementary excitation 
spectral density 

3 K 

= --lmGnu(E) = -llm == = = =(-l-22) 

I" Ti" vp - ft + *r) 2 - (J 2 + L> 2 ) 

It will be used in the next Section for examination of thermodynamical 
properties of the model in question. 

In order to estimate spin correlations one should calculate average 
fermion correlation function < Cm(Q)c n (t) > that can be found from the 
relation il 



i t_ r ^ i .•.^ Gf ™»( w +* c )/ 



< c+ (0)c n (i) > = Im / dwexp (-iwi) -^-r H.1.23) 

I" J-oo exp (Pw) + 1 

In particular, the static average fermion correlation function at low tem- 
pereture limit comes out from the following calculation 



< ct(0)c m+p (0) > = < cTnC m +p > = < C+_ p (0)c„(0) > 
1 f° — 

= Im / duG n (w + «e) 

71" J-oo 

7>" J-oo 



1 



71" |P 



Im i exp (itpp) y + ij- \/(y + ij) 2 - 1 |jj= 



1*1 .v=-n 



-oo 



TT \p\ 



Im {exp (itpp) 



+i"f — i 




V{uj 2 _ 7 2 _ 1)2 + 472w 2 + W 2 _ 7 2 _ J 



br 



, / V^l - 7 2 - I) 2 + 4 7 2 ^ 2 - o;g + 7 2 + 1 



,(1.24) 



where y = 



u>o, w 



— On 



f, 7 



= e+r 



The derived 



results (1.15), (1.21), (1-22), (1.24) remind the corresponding expressions 
obtained in slightly different cases in Refs. p4 Ea. 



2 Thermodynamical properties. The influ- 
ence of Dzyaloshinskii-Moriya interaction 

The obtained in the previous Section results are of great use in un- 
derstanding the thermodynamical properties of the model in question. 
Really, consider a model with a certain realization of transverse fields 
at sites. After exploiting Jordan-Wigner transformation (1.2) for the 
Hamiltonian (1.1) one comes to a quadratic in Fermi operators form that 
can be diagonalized by a canonical transformation r\ k — X)i=i(9fci c i + 

hkjCJ) § (see also (§, ||) with the result H = J2k=i A k(VkVk - §)• 
Elementary excitation spectrum A^ and the coefficients g k j and h k j are 



^JY 



N 



determined from A k g kn = J2i=i 9kiA in , -A fc /i fc „ = £\ =1 h ki A* in , where 



Aij = (fi + £U)Si; 



J+iD . 



J-iD 



energy per site for this realization is straightforward 



5j,,_i. The calculation of free 



/= lim — I --InTT 

= — - / dEp{E)\n [2 cosh 



exp 



/3A fc 



exp 



/3A S 






p(£) = lim l£«5(£_A fc ) 

iV— >-oo iV ^ — ' 



(2.1) 



and the result of its averaging over configuration is given by 
7 = - ~ J dE^E) In (2 cosh ^J 



= ~J dEp(E) In [1 + cxp {-(5E)\ - ^, (2.2) 



where p(E) is the average spectral density that has been found in the 
Section 1 (formula (1.22)). 

Further one finds the internal energy 

* = 1 + 4 = J dEm TT^m-T- < 2 ' 3 > 

the entropy 



* = %3 



J dEp ^ { ln t 1 + CXP H^l + l + exp(^) fc 4 > 



and the specific heat 



c = -/9-ii = .# / dEp(E)- -—-. (2.5) 



Using magical property of p(-E) (1-22) 



— ~TF)- — 



1 

-— Im- 



7T V(^-«o + ir) 2 -(j 2 + ^ 2 ) 



-^P(^) (2.6) 



one finds the transverse magnetization 

exp {(3E) 



<^h' > '^'! dEm T^W)-\ <17) 



3 = 1 

and static transverse susceptibility 



d< ± V s z > 



Note, that since at T — ► ti — '.aw ^ if £ > and t- — x ,„„, —> 1 if 

l+cxp {pE) 1+exp {pE) 

E < 0, the transverse magnetization and static transverse susceptibility 
at Qq = in low tempereture limit are given by 



1 1 

(2.9) 



< — > sf > 

N ^ J 



-Q 

dE' 
o 



1 
Im- 



7T v /( J B'+2r) 2 -(J 2 +D 2 ) 



XlI=---=JL===. (2.10) 

Ti" Vr 2 + j 2 + d 2 

The results of numerical investigation of thermodynamical proper- 
ties of the model (1.1) are depicted in Figs. 1-5. The dependence on 
Dzyaloshinskii-Moriya interaction (D — 0, D — 0.5 J, D = J) of the 
average spectral density p(E) (£! = 0) for different values of the width 
of lorentzian distribution T (T = 0, T = 0.25J, T = 0.5 J, r = J) is 
shown in Fig.l. In Figs. 2, 3 it is shown the temperature behaviour of en- 
tropy and specific heat for Slo = for different values of Dzyaloshinskii- 
Moriya interaction (D — 0, D — 0.5J, D = J) and different values 
of r (r = 0, T = J). The changes in the dependence of transverse 
magnetization on the value of transverse field tto that are caused by 
Dzyaloshinskii-Moriya interaction (D = 0, D = 0.5 J, D — J) for sev- 
eral values of T (T = 0, r = J) at T = can be seen in Fig. 4. In 
Fig. 5 the temperature behaviour of static transverse susceptibility in 
zero transverse field for few values of Dzyaloshinskii-Moriya interaction 
(D = 0, D = 0.5J, D = J) and V {Y = 0, T = J) a re pre sented. 

As it can be easily seen from the formula for p(E) (1-22) and in 
Fig.l the presence of Dzyaloshinskii-Moriya interaction formally causes 
the change of J 2 to J 2 + D 2 . This leads to effective increase of in- 
terspin interaction and results in broadening of zones both with sharp 
edges (when T = 0) and with smooth ones because of randomness (when 
1 ^ 0). Rather small quantitative changes in temperature behaviour 
of entropy and specific heat (Figs. 2, 3) reveal the tendency caused by 
Dzyaloshinskii-Moriya interaction. Dzyaloshinskii-Moriya interaction 
decreases the transverse magnetization for a given value of transverse 
field in both cases when T = and when r ^ (Fig. 4) and decreases 
the value of static transverse susceptibility in zero transverse field in low 
temperature region. The obtained result are in agreement with the de- 
rived ones for non-random version of II? s = \ isotropic XY model in 
transverse field (il| |l^] . 

Unfortunately, the obtained results do not permit to get any exact 
estimations for static spin correlation functions. Considering the sim- 
plest one < s Z jS z - +n >, using the relation s| = cj~Cj — | an d exploiting 
Wick-Bloch-dc Dominicis theorem one finds 



S 3 S 0+n > 



— + — + 1 — + — 1 — + 1 

= < CJC 3 cj +n C J+n > - -< CjCj > - -< Cj +n c j+n > + - 



< CJCJ >< C+ +n C j+n >-< C+C+ +n >< CjCj+n > 



P(E) 



0.4 




r = o 



0.2 



P(E) 



0.4 



2 ^ 



r = o.5J 



P(E) 



T = 0.25J 




0.4 



r = j 




Figure 1: The average spectral density p(E) (1.22) vs. ^ for different 
values of T; ilo = 0, D = (dashed curves), D = 0.5 J (long dashed 
curves) and D = J (solid curves). 
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r = o 




Figure 2: Temperature dependence of entropy s for T = and T = J: 
Oo = 0, D = (dashed curves), D = 0.5J (long dashed curves) and 
D = J (solid curves). 
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Figure 3: Temperature dependence of specific heat c for T = and 
r = J; fi = 0, D = (dashed curves), D = 0.5J (long dashed curves) 
and D = J (solid curves). 
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Figure 4: The dependence of transverse magnetization m z = 

< W EJLi s ] > ° n transverse field ^ for T = and T = J; 1//3 = 0, 
Z? = (dashed curves), Z? = 0.5J (long dashed curves) and D = J (solid 
curves) . 
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Figure 5: Temperature dependence of static transverse susceptibility \zz 
for T = and T = J; Q = 0, D = (dashed curves), D = 0.5 J (long 
dashed curves) and D — J (solid curves). 
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-< cfc j+n >< c jC + +n >-< c+c, > + -. (2.11) 



Thus, knowing only < CmC n > it is not possible to derive any rigorous 
results even for simplest (transverse) static spin correlation function. It 
is worthwhile to note that its calculation faces with the similar problems 
that were considered in Jl| . Apparently, the problem of spin correlations 
can be solved within developed recently numerical approach f2q , [27| , |28| , 

HE§. 

3 Standard approximate approaches in the 
theory of disordered spin systems 

Let's consider the results one faces with after adopting some standard 
approximate approaches in disordered spin systems theory while consid- 
ering the thermodynamics of the model defined by (1.1). 

3.1 Bose commutation rules approximation for spin 
operators s + , s~ 

Implying instead of Pauli commutation rules for spin operators s + , s~ 

S 7> 5 ml =S jm (l~2.S+SJ 



i s j ' s m] = <5jm(l — ^ s t s j ) Bose commutation rules, that is 



[Sj.il-^™, (3-1) 

one comes to a system of bosons on lattice that may transfer from site to 
site with random (lorentzian) energy at sites. Introducing the following- 
Green's functions 

D* m (t) = TiO{±t) < [a-(t),a+(Q)] >, (3.2) 

and repeating the derivation of Section 1 one ends up with 



D^ m (u±ie) 



f- u-Q ±i(e+r) _ /r ^-On±i(e+r) -i 2 ~~n \n-m\ 

V J 2 + D 2 / r ^-n ±i(6+r) 12 , 



VJ 2 +D 2 

that gives for average spectral density formula (1-22). However, the 
average boson correlation function < Sm(0)sn(t) > = — — ImJ"_ diox 
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exp (— iu>t) 



D nm (uj+ie) 
cxp (/3cj) — 1 



contains the denominator that tends to when 



Let's discuss this problem in details. After assuming Bosc commuta- 
tion relation one has the quadratic in Bose operators s'f, sj form (1.1) 



iY 



that with the help of linear canonical transformation 7^ = E 7=1 (/< 



I-]' 



dkjst) can be diagonalized with the result H = — \ ^7=1(^0 + ^ 



,A r 



Sfe=i ^kjtlk + const. Since £k, fkj, dkj are determined from the equa- 
tions Skfkn = Ei=i /fei^m, ~£kdkn = J2i=i d ki A* n , one immediately 
concludes that £t — Afc. The ground state energy of this Hamiltonian 



coincides with the exact value if const = — i X^fe=i ^-k 



^N 



iEjLiCV 



n. 



and thus finally i7 = X)fe=i ^fe(7fc 7fe 
given by 



2). Free energy per particle is 



f = lim — x 

JV^oo AT 

J --In J] exp 



all At > 



lim 

JV->oc 



/3A fc 



1 

N 



[1 + cxp (-/?A fe ) + cxp (-2/?A A 



5>T 



exp 



2 



- exp (-/3A fe ) 
dEp(£) In [1 - exp (-/?£)] - - / dEp(E)E 



(3.4) 



Note, that the partition function for Bose system exists only when the 
condition all A K > (or p(E) = if E < 0) is valid. Really, otherwise 
the state without bosons is not the ground state, since one-boson state 
has smaller energy, the energy of two-bosons state is more smaller etc., 
the probability of their appearance increases respectively and thus the 
partition function tends to infinity. This difficulty arises because of ap- 
proximate treating of elementary excitations as bosons (they are exactly 
fcrmionic objects in the case under consideration) and it is crushing if 
p(E) ^ for E < 0. Thus, one can consider the average thermodynami- 
cal quantities obtained within approximation (3.1), that is, free energy 



/ = 



dEp(E) In [1 - cxp (-/3E)} - 



n 



internal energy 



dEp{E) 



E 
exp (J3E) - 1 



2 ' 



(3.5) 



(3.6) 
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entropy 

(3E 



= J dEp(E) | - In [1 - cxp {-[IE)] 



exp (fiE) - 1 
specific heat 

E 2 



(3.7) 



T ^ ldEp{E) (^^W (3 ' 8) 



transverse magnetization 



<h^>=J dE ^i^m^- 1 * (3 - 9) 



and static transverse susceptibility 



-- = - s, ! dmE] i^km (3 - 10) 



(2sinh^) 
when the relation 



p{E) = for E < (3.11) 



holds true. 

Considering at first the case T = 0, when 



p{E)={ ^■P+D--(E-n a )^ v v ( 3 . 12 ) 

I 0, otherwise, 

one finds that (3.11) is valid only in the case of strong transverse fields 



ft > VJ 2 + D 2 ; (3.13) 

Dzyaloshinskii-Moriya interaction increases the value of the field behind 
which (3.11) becomes true. In the case T ^ one immediately finds that 
(3.11) is never true, since there always will be elementary excitations 
with negative energy and therefore it is impossible to treat them as 
bosons. Thus, the consideration of the disordered version of the model 
demands the revision of the problem of validity of approximation that 
was suitable for study of non-random version of the model. 

The results of numerical calculation of temperature behaviours of 
entropy and specific heat according to exact formulae (2.4), (2.5) (solid 
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lines) and approximate ones (3.7), (3.8) (long dashed curves) in the case 
of validity of approximation (3.1) (r = 0, Oo = (\/2 + 0.1) J, D = 
0, D = 0.5J, D = J) are presented in Figs. 6, 7. These results show 
that Bose commutation rules approximation for spin operators s + , s - 
gives suitable results only for low temperatures, and in the presence of 
Dzyaloshinskii-Moriya interaction only at very low temperatures. 

3.2 Tyablikov-like approximation 

Does not assuming Bosc commutation rules for operators s + , s~ (3-1), 
one faces with the equations of motion that contain more complicated 
Green's functions 



ij t DZ n (t) = 6(t)6 nm (-2 < < >) + (fi + n n )Dt m (t) 



J + iD 



(-2)( T iO(±t)<[s z n (t) S - 1 (t), S +}>) 



2 
+ ^R(-2)( T id(±t) < [s-_ 1 (t)sUt),s+] >)• 



(3.14) 




Figure 6: Entropy vs. temperature for r = 0, Qq = {\/~2 + 0.1) J: ex- 
act results (solid curves), Bose commutation rules approximation (long 
dashed curves), Tyablikov-like approximation (dashed curves). 
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Within Tyablikov-like approximation it is supposed that 

< s z n (t)s- +1 (t)s+ >~ <^^ < s- +1 (t)s+ >, 

< s+s z n (t)s- +1 {t) >« <^^ < s+s- +1 (i) >, 

< s-_ 1 (t)s z n (t)s+ >« <s z > < s-_ 1 (t)s+ >, 

< *m«n-l(*)*n(*) >~ < S " > < s lVlW > ■ (3-15) 

Then instead of (3.14) one has 

ij t Dl m {t) = -2<^>S(t)5 nm + (fio + nn)I>^ m (t) 
+(-2<^>)^±!^^ +lim (t) + (-2<^>)^^^_ 1;ro (t).(3.16) 
Acting like in Section 1 one ends up with 



D% m (u)±ie) = 

r coi-n ±i(e+r) /r M-S"i ±i(6+r) 12_-|\I"- TO I 

exp [iy>(ra— m)] i (-2<s*>)VJ 2 +-D 2 V [ (-2<s">)V./ 2 +-D 2 J * 

VJ 2 +£» 2 /r "-"o±»(*+r) 12 1 

V L (-2<^>)V./ 2 +c 2j 

that yields the following result for average spectral density 

1 1 



(3.17) 



p{E) = — Im == — (3.18) 

The introduced average transverse magnetization at site is determined 
from the equation 



< s z > = < s+s > = Im / duo 



1 1 T Z" 00 , D mm (uj + ie) 1 



2 7T ./_„, exp(/?w)-l 2 



, , - p ( w )_-. (3.19) 

L^exp^)-^ 1 ; 2 

Equation (3.19) contains Bose factor ,* , 1 and thus, apparently 

Tyablikov-like approximation is possible if p(u>) = for uj < 0. The tem- 
perature behaviours of entropy and specific heat found within Tyablikov- 
like approximation (3.15) (dashed lines) in comparison with exact re- 
sults (solid lines) and results obtained within approximation (3.1) (long- 
dashed lines) are shown in Figs. 6, 7. Despite some improvmcnt over the 
Bose commutation rules approximation that can be seen in Figs. 6, 7 the 
results obtained within Tyablikov-like approximation generally speaking 
are not closer to exact ones in comparison with the results derived within 
Bose commutation rules approximation. 
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3.3 Randomly disordered crystals theory methods 

Let's consider ID spin-^ isotropic XY model with Dzyaloshinskii-Moriya 
interaction in random (not necessary lorentzian) transverse field within 
usually used methods in the theory of disordered crystals || |], [51j . The 
starting point is the equation for Green's functions that after introducing 
the notations 

W pr = £lp8 pr , Vp S = <$;s,p+l H n Ss : p-1 (3.20) 

can be written in the form 

wG£ m (w) = S nm + tt Q G* m (uj) + W nr Gf m (co) + V ns Gj m (u®.21) 

(for Gj m (w) defined by (1.4) or (3.2) within Bose commutation rules 
approximation) or 

uD* m {u) = -2 < < > S nm + n D* m {u) + W nr Df m (Lo) 

+ (-2 < < >)V ns Dj m (u;) (3.22) 

(for D^ m (uj) (3.2) within Tyablikov approximation, < s z n > is only ther- 
modynamically averaged value (without configurational averaging as in 
(3.15)) of transverse spin at site n); the summation over the repeating 
indices from 1 to N is implied. The different approaches of randomly 
disordered crystals theory are constructed from so called propagator and 
locator expansions. 

Propagator expansion. Let's rewrite the Hamiltonian of the system 
in question (1.1) in the form 

N N , . 

h = a H + J2 a,*! = ci h + E n A 4 s J 2 (3 - 23) 

3=1 3=1 V / 

and introduce Green's functions a G^ m (cu) or °Djf m (w) for the system 
with Hamiltonian a H that, naturally satisfy the following equations 

u a GT m (uj) = 5 nm + n a Gf im {u;) + V ns a Gj m (u;) (3.24) 



or 



u a D* m (u>) = -2 < a* > S nm + n a D* m (u>) 

+ (-2 < < >)V ns a Dj m (u). (3.25) 
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Multiplying (3.21) (or (3.22)) by a GJ n {cj) ( a DJ n (uj)) one comes to the 
following equations 

G*„» = a G* m {u) + a G* (w)W nr G£» (3.26) 

or 

D£»= a J Dj m H^|^+ -I>T(u;)W' nr I>T i (a,), 

Wnr = Wnr z + V nr ( ^4^ - ^ • ( 3 - 27 ) 

-2 < s z > \<s z > J 

Note, that since generally speaking < s z n >^< s z > even in the case of 
diagonal disorder after Tyablikov-likc approximation one faces with non- 
diagonal disorder problem. Expanding of (3.26) (or (3.27)) in degrees of 
W nr (W n r) leads to propagator expansion. 

Locator expansion. Let's introduce locators 

ffnm( W ) = 9n S nm = /o , n \ 6nm ( 3 " 28 ) 

U) — [}Iq + \l n ) 

or 

<£ m (w) = dl8 nm = ~ 2 * **"„ J nm . (3.29) 

Then equations (3.21), (3.22) can be rewritten in the form 

Gnm(v) = 9nm(u) + 9np(^WpsGj m (uj) (3.30) 

or 

D* m (u) = d&Jw) + dl p {Lo)V ps Dl n {Lo). (3.31) 

While expanding r.h.s. of (3.30), (3.31) in degrees of V ps one comes to 
locator expansion. 

Further analysis deals only with diagonal disorder when, for example, 
the propagator expansions have the form 

G* m (u>) = a Gj m (co) + 
a Gf n (cj)n n a G* m (w) + a G* n (w)n n a GT p (u,)n p a G^ m (cj) + ... (3.32) 

or 

DIM = °l£» + 
a Df n (aj)Cl n a D* m (u;) + a DjJco)n n a D^ p (co)n p a D* m {u) + ... , 

fi n = -Cl n /2< s z >.(3.33) 
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Extracting in (3.32) ^-matrix 



a 



one can rewrite (3.32) as a series in degrees of i-matrix 

GlnH = a GT ( W ) + -GT MO, a G* m (. 



(3.34) 



+ G£»ft„ a G*»ft„ a G* m (u; 

+ a Gf n (co)n n a Gl p (n ^ p) {io)% »Gj m ( W 

+ °G^(a;)a, G^n(w)0, "G^MO, a G^ ? > 

+ °G* (w)fi n Q G^ p {n ^ p) (co)a p G£,(w)fi p a GT m {u) 

+ °G* MO, a G* n (u)Cl n a Gl f (n#/) (u;)fi, a Gj m (co) 

+ a G*»o„ a ci p {n ^)% a c; f (p#/) H% °g£» + ... 

= »Gj ro ( W ) + °G*»i„ a G* ro (u;) 
+ a G* n (u)t n a Gl p {n ^ p) (io)t p a GT m (co) + ... . (3.35) 

Within approximation of average i-matrix one assumes that 

- -a 



t n — t = / d£li...dQ,Np(£li, ••■, &n) 



1- a G%„(u>)flr, 



(3.36) 



where 



a G* m (u) = 

exp [itp(n — m)\ ^ ^.p+d 2 



//■ ui — Qg±ie \2 illw— ml 
V ^ V.7 2 +-D 2 ' J 



V J 2 + £ 2 



V^ 



a; — Op±2e \2 



- 1 



(3.37) 



VJ 2 +£> 2 
In result one is able to sum the series for average Green's functions (3.35) 

°G T H \ 
.1- «G*{w)t) gm 



Gj m (u) = 



(3.38) 



Within coherent potential approximation one should seek the Green's 
functions of the system with Hamiltonian 

H = 

N JV-1 JV-1 

£(fi + fya'j+J Y, (***j+i+* v j* v J+ i)+D E (s** V j + i-* V j*i+i)> ( 3 - 39 ) 
J— 1 J— 1 J =1 
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where Cl is unknown coherent field. These Green's functions are given 

by 



t u — Qg — n±ie _ I ( u — Op — i1±ie ^2 _ il|n— m| 

cxp [iip(n - m)\ 1 772qrp2 y I y, 7 2 +g 2 ; -H 

V J 2 + D 2 /( u-n -n±ie \ 2 -, 



Since 



N ( 1\ 

H = H + Y / (n j -(l)lsp---), (3.41) 

j=i \ / 

w(^ m («) - <f nm + (fio + A)(^ m (w) + K S G^», (3.42) 

one can get acting like while deriving (3.26) the following equation 

Gf m (") = Gf rn (uj) + G* n {w)W nr G? m {u), 

W nr = W nr - (I5 nr = (fl n - Ci)S nr , (3-43) 

and hence the following propagator expansion 

Gj m (u;) = Gf m (cj) + (3* (w)(0, - ft)(5* m (w) 

+£* (w)(fi n ~ n)GT p (Lo)(n p - ft)G*» + ... . (3.44) 

This series can be rewritten as an expansion in degrees of t -matrix 

i ll n ~ \L . . 

"~ l-Gf n (u;)(fl n -ny 

namely, 

G$» = 
^M+^M^^ m M+^(w)tn^ p(n ^ ) M* P <^ n M+...(3.46) 

Determining the coherent field O from the condition 

T n = [dn 1 ...dn NP (n 1 ,...,n N ) """" — -r = 0(3.47) 

where accorging to (3.40) 

Gl n (w ± te) = 1 (3.48) 

/ (w-fi -^±«e) 2 -( t / 2 +£' 2 ) 



21 



one finds that 



Gf m (w) = <2£» + 6*(u)i n Gl m {w) 



+G* n (w)i n G* p {n ^ p) (cj)i p GJ m (cj) + ... ~ G* m (uj), (3.49) 

that is the desired result within coherent potential approximation. 

In the case of lorentzian transverse field the equation for coherent 
field il (3.47) reads 

1 r 

dfL — - 



7r(fij +?r)(fij--ir) 

(£L - f2) A /(w - fi - ^±«e) 2 - (J 2 + £> 2 ) 
/ ; V = = = 0. (3.50) 

(w - fi - fi±ie) 2 - (</ 2 + £> 2 ) - Oj + & 



Supposing that Im[£l + J (lo — Oo — f2±ie) 2 — (J 2 + D 2 )] one can 

perform integration with the help of the residuum theory getting in result 
instead of (3.50) 



(=F»r - fiW(w - fio - 0±ie) 2 - (J 2 + D 2 ) 

= 0. (3.51) 

(w - fi - &±ie) 2 - (J 2 + D 2 ) + n±ir 

Equation (3.51) has solutions Cl = =pir and after insertion them into 
(3.40) one gets exact result (1.21). It can be proved post priory the pos- 
sibility of assumed displacement of poles in (3.50) at least for u> — > oo. 
Thus, in the case of lorentzian transverse field the coherent potetial ap- 
proximation contains exact result for average Green's functions G| m (w). 
Consider now another version of random transverse field that is given 
by probability density 

N 

p(fii, ...,n N ) = Y[[x6(ilj) + (1 - x)6{ilj - 13)}. (3.52) 

Then the equation for coherent field Cl (3.47) reads 

~ L? (1-x) J^ ^=0 (3.53) 



i - GT n (uj)(-n) ' i - g*„(u])(v - n) 
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and after some calculation reduces to 3th order algebraic equation for tt 



n 3 



{x 2 -2x)l5 2 -(J 2 +D 2 )+4(l-x)l5{u} -n ) + (uj-n ) 2 ^ 2 



2(xl5+tt Q -Lu) 
(l-x)l5(J 2 +D 2 )-(l-x) 2 l5 2 (uj-n )-(l-x)l5(uj-n ) 2 



a 



xl5+flo— lu 

-(1 - x) 2 l3 2 (J 2 + D 2 ) + (1 - x) 2 U 2 {oj - n ) 2 
2(xV + fl - to) 

= Q 3 + a(l 2 + bQ + c = 0. (3.54) 

It is generally-known [53 that at first one should substitute Cl = y — | 

obtaining in result y 3 +py + q = with p = — ^- + b, </ = 2 (|) — ^ + c. 
Then the real " noncomlete" cubic equation has 

• one real and two conjugate complex roots if Q > 0; 

• three real roots at least two of which coincide if Q = 0; 



• 



three different real roots if Q < 0; 



here Q = (§) + (|) . The solution can be presented in the trigonomet- 
rical form 

• if Q > 0, p > 0, then 

/P IP ( r \ 

-cot 2a, 2/2,3 = \ r I cot2a±iv3csc2a) , 



• if Q > 0, p < 0, then 



tana = ^/tan| (|o| < J) , tan/3 = -^(J) (|/3| < |) (3.55) 



i/i = — 2W — -esc 2a, J/2,3 = y — r (csc2a±«v / 3cot 2a j , 

y O y O \ / 

tana = ^J (|a|<£) , s inP=-J(-^) 3 (|/3|<|) ; (3.56) 

if Q < (and hence p < ), then 

„ / 23 a / p" /a 7r\ 

J / 1 =2^--cos-, y2 , 3 = -2^--cos(-±-j, 

cosa= q (3.57) 
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(all cubic roots are real). 

Further one should act in a following way: 1) for a given x to calculate 
a, b, c in (3.54), p, q and Q; 2) comparing Q and p with the zero, to 
write down 2/1,2,3; and hence fii.2,3; 3) to check whether t after inserting 
Clj is really equal to zero (since the algebraic transformation from (3.53) 
to (3.54) may lead to appearance of extra roots); 4) to substitute flj 
into G~ n (E + ie) (3.40); and 5) to calculate the average spectral density 
within coherent potential approximation (3.47) 

1, A ,„s 1, 1 



p(E) ~ --ImG- n (E) = --Im ; = . (3.58) 

V 



" l(E-no-n 3 +ie) 2 -(J 2 +D 2 ) 



Further the case J = 1, D = 0, = 1 is under consideration. Here 
the example of performing of this program is given: 



a; = 0.01, 

£ = 0.02, 

a = 47.015001, b = -97.000202, c = 48.985399, 
p = -833.803646, q = 9267.099096, Q = 0.896875, 
fii = -49.014405 + i0.000000, 

2 = 0.999702 + i0.001967, 
(l 3 = 0.999702 - i0.001967, 

|t|j 2i =1.012523xl0- 13 , |t|^ 2 =3.978199xl0- 3 , |t|^ 3 =3.064551 x 10" 

E = 0.03, 

a = 22.505001, b = -47.985301, c = 24.480448, 

p = -216.810318, q = 1228.762879, Q = 0.087101, 

Oi = -24.504032 + i0.000000, 

n 2 = 0.999516 + i0.002358, 

q 3 = 0.999516 - i0.002358, 

|t| fii = 1.187939xl0- 14 , |t| ft2 =4.811497xl0- 3 , 1*1^=2.295226x10" 

E = 0.04, 

a = 14.331667, b = -31.640401, c = 16.308864, 

p = -100.105961, q = 385.512229, Q = 0.023724, 

S7i = -16.330344 + i0.000000, 

(l 2 = 0.999339 + i0.002665, 

3 = 0.999339 -i0.002665, 
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|t| r - 2i =5.717649xl(r 15 , |% 3 =5.488112xl(r 3 , |i|^ 3 =2.070775xl0- 12 ; 

£ = 0.05, 

a = 10.242501, b = -23.463002, c = 12.220623, 

p = -58.432610, q = 171.921947, Q = 0.009691, 

fii = -12.240840 -M0.000000, 

2 = 0.999170 -M0.002918, 

ft 3 = 0.999170-^0.002918, 

|t| fii =1.362799xl0- 14 , |% 3 =6.062763xl0- 3 , |£|^ 3 =3.098188xl0- 12 , 

etc. Since t with Cl\ for x = and x = 1 is not equal to zero, this 
root has been rejected. Unfortunately, there are no exact results for 
a random model in question (3.52). Nevertheless, for arbitrary random 
spin-1 anisotropic XY model in transverse field it is possible to calculate 
numerically (see §|, |7)) the quantity R(E 2 ) = ± £f=i S(E 2 -A 2 k ), the 
average value of which is connected with the average value of p(E) by 
the relation 



p{E)+p{-E) 



W&) = 2| - '- 0.59) 

In Fig. 8 the results of calculations within coherent potential approxi- 
mation (3.58), (3.59) (broken lines) are depicted together with the exact 
results (solid lines). The comparison shows just how little is the change 
in R(E 2 ) and thus in thermodynamical quantities. This seems to be 
conditioned by the fact that for the model in question (since it is de- 
scribed by Hamiltonian (1.3)) the thermodynamical averaging has been 
performed exactly. 

4 One-dimensional spin-| XX Z Heisenberg 
model with Dzyaloshinskii-Moriya inter- 
action in lorentzian random external field 

This Section is devoted to examining of somewhat more complicated than 
(1.1) case. Namely, now the intersite interaction contains the coupling 
of z components of neighbouring spins and the Hamiltonian reads 

N JV-1 JV-1 
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N-l N 

+ JZ E *i*i+i = E("o + ^) 

N-l 



4 S J 2 



E 

3 = 1 



J + iD 



• +o- 



S 7 S 7 + l + 



J-iD _ + 

— 2~- S 3 S 3+l 



N-l 
3 = 1 



S 3 S 3 S j + l S 3 + l 2 S 3 S 3 2 Sj+1 ' Sj+1 



I M 4 - 1 ) 



After Jordan- Wigner transformation (1.2) one gets 



H = - 



NQ (N-1)J z 1 



JY 



-3l> 

3 = 1 



J z 



Qo + Qj - ^— (2 - 5j,i - 5 j>N ) 



N r 

E 

^ fJ + iD , J-iD 

M - r - c i c i+i — — 



ci Cj + 



-j^i+i 



N-l 



d'l^i ) +j z 2J c -' c ' r 



j" < -J l -.7+ 



iCj+{4.2) 



3 = 1 



Note, that the key difficulty is that due to new intersite interaction one 
faces with the terms that are the products of four Fermi operators. It is 
easy to find the equations of motion for Green's functions (1.4) 



ijGl m {t) = 5{t)5 nm + (n +n„- J z )Gl m {t) 
J + i-D^- J ' — iD x 

+ — 2~ G n+l,m(*) + — 2~ G n-l,m(*) 

+ J*[T*«(±*) < {( CnC + +1 C„ +1 )(i),C+(0)} > 

T»0(±*) < {(c+_ lC „_ lC „)(i),c+(0)} >]. 



(4.3) 



Equations (4.3) contain higher Green's functions and thus cannot be 
solved exactly. Nevertheless, making the approximation 



< {( C n+l C n+lCn)(t),C^(0)} >RJ < C+ Cj > < {c n {t) , C+} > 
_< C 3 + 1 C 3 > < { C ™+l(*)' C m} >, 



< {(c+_ 1 c„_ 1 c„)(i),c+(0)} >RJ < C+ Cj > < {c„(i),c+} > 



■<CJ_ lCj ><{c n -l(t),C+}>, 



(4.4) 
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one gets instead of (4.3) the equations 



ijGl m (t) = 5{t)5 nm + [Sl + J z (2< c+c 3 > - 1) + Sl n ]G? im (t) + 



J-2J Z < c+ Cj >+iD T J-2J Z < cj ,cj >-iD T 

o u n+l.ml l i+ o "n-l,'. 



.(*)(4.5) 



that can be treated like (1.6). Substituting SIq + J z (2< cJcj > — 1) 

instead of Slo, J — 2J Z < ct ±1 Cj >±iD instead of J ±iD into (1.15) one 
finds that 

1 



GT(uj±ie) 



Sl + J+ CX P (* K ) + J- exp {—in) ±i(e + T) 



Slo = Sl + J z (2<cJ Cj >-l 



J+ = 



J - 2J Z < ct ±1 Cj >±i£) 



(4.6) 



Returning back to the site representation needs the calculation of the 
following integral (compare with (1.16)-(1.21)) 



i r 

Gn m (w±ie) = — / eta exp [i(n — m)n\ 
2lT J-k 

1 



Slo + J+ exp (in) + J- exp (— in) ±i(e + T) 



(4.7) 



Putting z = exp (±in) (for n > m and n < m, respectively) one comes 
to contour integral over unit circle centred at z — in complex plane z 

1 



G^ m {oj±ie) = --— f dz = 

2m J J+z 2 - [w - Slo±i{e + T)]z + J_ 

dz- 



2m J J+(z — Zi)(z — Zj ) 



U) - 


~ Sl ±i(e + T) + yj[u - £7 ±z(e + T)] 2 - 


- 4J+ J_ 




2J+ 




u> - 


- I^±i(e + T) - \J[uj - H^±i(e + T)] 2 - 


- 4J+ J_ 



2J, 



Z l Z 2 



+ J_ J - 2J Z < cJ^Cj >-iD 
J + J - 2J Z < cf +1 Cj >+iD 



(4.8) 
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for n > m, 



Gtm{uJ±ie) 



1 
"2kI 



dz- 



y \n— m\ 



J-z 2 -[u- fl ±i(e + T)]z + J+ 



1 

2^i 



dz- 



-\n—m\ 



J-(z-z 1 )(z- Z 2 ) 



U! - 


- n ±i(e + r) + \J[u> - n ±i{e + r)] 2 - 


-4J+J_ 




2J_ 




UJ - 


- H^±j(e + r) - yj[u - W Q ±i(e + T)} 2 - 


-4J+J_ 



2J_ 



J -2J Z < cj +1 Cj > + iD 



J- 



J - 2J Z < c+_ x Cj >-iD 



(4.9) 



for n < m. The result of integration can be easily obtained; it depends 
on the position of z 1 2 in complex plane z with respect to unit circle 
centred at the origin of complex plane. The introduced static average 
fcrmion correlation functions are given by 



< C+Cj > 



< Cj-lCj > 



1. 



-Im 



du> 



GjAu + ie) 



exp (/3u>) + 1 



(4.10) 



< 4 C 3 + 1 > 



1 

-— Im 

7T 
1. 



duJ ln \ 

_ x exp (fiuj) + 1 



(4.11) 



< c t+i c i > = < CjCj-x > = --Im 



dw G Ii±^ + i ^ 



(4.12) 



71" J-oo ex P (A*>) + 1 

The equations (4.8)-(4.12) permit to determine self-consistently cor- 
relation functions < cJcj >, < c] +1 Cj >, < c^-^Cj > that were intro- 
duced in (4.4). Knowing them one can get the average spectral density 
that is given by 



p(E) = 



— Im — 

7T 2-7T 



dn- 



ilo+Jjt- exp (in)+ J_ exp (—in) +i(e+F) 



— Im 

7r 2ni 



dz 



T- ( 413 ) 

J±(z-Zi ){z- z 2 ) 

Such approach can be viewed as a generalization of rather old ideas (see 
[p3| ) for the case of random lorentzian external field. 
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5 Conclusions 

The results which have been obtained in the present paper are summa- 
rized as follows. Considering spin-i isotropic XY chain with Dzyaloshin- 
skii-Moriya interaction in random lorentzian transverse field (1.1) as a 
system of fermions (1.3) one is able to obtain exactly the average one- 
fcrmion Green's functions G„ m (w±ie) (1.21). This permits to get exact 
result for average spectral density (1-22) and therefore, for thermody- 
namical quantities (2.2)-(2.5), (2.7)-(2.10), to examine the influence of 
Dzyaloshinskii-Moriya interaction on them (Figs. 1-5). Unfortunately, it 
appeared impossible to reach any exact results concerning spin correla- 
tions in model in question. The system that is considered does not rep- 
resent accurately any physical system, but may serve as a testing ground 
for theory and numerical techniques and in this sence, has proved most 
valuable. The results of comparison of conclusions obtained within differ- 
ent approximate approaches with exact ones (Figs. 6, 7, Fig. 8) represent 
some interest and may be useful in understanding the region of validity 
of various well-known methods of randomly disordered crystals theory. 
After introducing into the Hamiltonian of the interaction between z com- 
ponents of neighbouring spins one is unable to obtain exact results since 
these terms lead to appearance of products of four Fermi operators in 
(4.2). After adopting approximation (4.4) the problem reduces to al- 
ready considered one, however, more complicated because of necessity 
to determine self-consistently three fermion correlation functions. 

The performed exact calculations for the system with Hamiltonian 
(1.1) or (1.3) are anticipated to be useful in leading to an understand- 
ing of the properties of one-dimensional disordered models of condensed 
matter physics. 
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Figure 7: Specific heat vs. temperature for T = 0, f^o = (V2 + 0.1) J: 
exact results (solid curves) , Bose commutation rules approximation (long 
dashed curves), Tyablikov-like approximation (dashed curves). 
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Figure 8: R(E 2 ) vs. E 2 : exact results (solid curves) (the averaging is 
done only over few random realizations) and the results within coherent 
potential approximation (long dashed curves) for the model with disorder 
(3.52): O = 0, J = 1, D = 0, = 1. 
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